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ABSTRACT 

The standard genetic code is used by most living 
organisms, yet deviations have been observed in 
many genomes, suggesting that the genetic code 
has been evolving. In certain yeast mitochondria, 
CUN codons are reassigned from leucine to threo- 
nine, which requires an unusual tRNA Thr with an 
enlarged 8-nt anticodon loop (tRNA| hr ). To trace its 
evolutionary origin we performed a comprehensive 
phylogenetic analysis which revealed that tRNA{ hr 
evolved from yeast mitochondrial tRNA H,s . To 
understand this tRNA identity change, we per- 
formed mutational and biochemical experiments. 
We show that Saccharomyces cerevisiae mitochon- 
drial threonyl-tRNA synthetase (MST1) could attach 
threonine to both tRNA| hr and the regular tRNA^ hr , 
but not to the wild-type tRNA His . A loss of the first 
nucleotide (G_ 1 ) in tRNA Hls converts it to a substrate 
for MST1 with a K m value (0.7 \iM) comparable to that 
of tRNA{ hr (0.3 nM), and addition of G_! to tRNA| hr 
allows efficient histidylation by histidyl-tRNA 
synthetase. We also show that MST1 from Candida 
albicans, a yeast in which CUN codons remain 
assigned to leucine, could not threonylate tRNAl hr , 
suggesting that MST1 has coevolved with tRNA{ hr . 
Our work provides the first clear example of a recent 
recoding event caused by alloacceptor tRNA gene 
recruitment. 

INTRODUCTION 

Transfer RNAs (tRNAs) are adaptor molecules that pair 
each amino acid with corresponding codons on the 
mRNA (1). Aminoacyl-tRNA synthetases (aaRSs) 



attach amino acids to the 3'-terminus of tRNAs (2), and 
the resulting aminoacyl-tRNAs (aa-tRNAs) are delivered 
by elongation factors to the ribosome, where codon-anti- 
codon recognition defines the genetic code of life (3). The 
genetic code was once thought to be universal and frozen 
(4), but later studies revealed multiple codon reassignment 
events in nuclear and mitochondrial genomes (5-7). These 
reassignments include stop-to-sense, sense-to-sense and 
sense-to-stop codon changes; they have been found in 
nuclear genomes from all three domains of life. In some 
bacterial, archaeal and eukaryotic species the UGA codon 
is reassigned to allow insertion of selenocysteine (Sec), the 
21st natural amino acid, in the presence of tRNA ec with a 
UCA anticodon (6,8). In a number of methanogenic 
archaea a UAG stop codon encodes the 22nd amino 
acid, pyrrolysine (Pyl), implemented by a tRNA Pyl which 
recognizes UAG (9,10). A major reassignment in mitocho- 
ndria is recoding UGA to tryptophan (Trp) implemented 
by a tRNA Trp with a mutated anticodon (7). Furthermore, 
the mitochondria of several yeast species reassigned the 
AUA codon from isoleucine to methionine (Met) by the 
abnormal recognition of AUA by tRNA Met (7). 

One last codon reassignment in the well-studied 
organism Saccharomyces cerevisiae is still not understood. 
In certain budding yeasts (including Saccharomyces, 
Nakaseomyces and Vanderwaltozyma) the four CUN (N 
denotes U, C, A or G) codons in the mitochondria are 
reassigned from leucine (Leu) to threonine (Thr) (11,12). 
This results from the loss of tRNA^ G (with a UAG anti- 
codon) that would translate CUN codons, and from the 
presence of an abnormal tRNA™ r with an enlarged 8-nt 
anticodon loop and a UAG anticodon (Figure 1) (11,13). 
This codon sense change was confirmed by protein 
sequencing of the S. cerevisiae mitochondrial ATPase 
(14). Mass spectrometry studies have also validated that 
at least three CUU and two CUA codons in S. cerevisiae 
mitochondrial-encoded proteins are recoded as Thr (15). 
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In addition to tRNA[ 11 , yeast mitochondria also express a 
normal tRNAj hr with a UGU anticodon that reads the 
ACN threonine codon box (Figure 1). It was shown pre- 
viously that while mitochondrial extracts from a wild-type 
(WT) 5". cerevisiae strain could attach Thr to both 
tRNA™ r and tRNAj hr , extracts from an MST1 mutant 
strain could only threonylate tRNAj hr but not tRNA™ 1 ' 
(16). Thus it was suggested that MST1 serves as a mito- 
chondrial threonyl-tRNA synthetase (ThrRS) specific for 
aminoacylation of tRNA[ hr , while a different ThrRS 
aminoacylates tRNA™ 1 in the mitochondria. Three 
decades after this discovery it is still a mystery how the 
unusual tRNAf hr emerged in the mitochondrial genome. 
Previous hypotheses suggest that tRNA}' 11 might have 
evolved from tRNAj hr , or alternatively from the missing 
tRNA[}™ G (11,17,18). However, both hypotheses lack 
convincing experimental evidence. To our surprise, bio- 
chemical and phylogenetic analyses demonstrate that 
tRNAf ir directly evolved from mitochondrial tRNA His . 
Saccharomyces cerevisiae mitochondrial tRNA Hls shares 
high (72%) sequence identity with tRNA™ 1 ', and a 
single-nucleotide change converts tRNA Hls to a substrate 
for MST1. Our work thus resolves the long-standing 
question regarding the origin of tRNA]" hr and its coding 
response. 

MATERIALS AND METHODS 

Cloning, mutagenesis and general methods 

Saccharomyces cerevisiae MST1, S. cerevisiae HTS1, 
Candida albicans MST1 and Schizosaccharomyces pombe 
MST1 genes were cloned into pET28a expression vector 
(Novagen) with an N-terminal six-His tag. Expression of 
recombinant proteins was induced at 37°C for 4h with 
0.5 mM isopropyl (3-D-l-thiogalactopyranoside (IPTG) in 
Escherichia coli strain BL21 -codon plus in Luria-Bertani 
(LB) media. His-tagged proteins were purified according 
to standard procedures. Mitochondrial tRNA genes were 



cloned into pUC18 vector (GenScript), and mutations 
were introduced using QuikChange Site-Directed 
Mutagenesis Kit (Stratagene). 

In vitro assays 

In vitro tRNA transcripts were obtained using the T7 
RNA polymerase runoff procedure as described (19). 
Aminoacylation experiments were performed as described 
(20) in the presence of 100 mM Na-HEPES pH 7.2, 30 mM 
KC1, 10 mM MgCl 2 , 2mM ATP, 25 uM [ 3 H] Thr (100 uCi/ 
ml), 25 uM [ 3 H] His (150uCi/ml) or 50 uM [ 3 H] Leu 
(lOOuCi/ml), 0.2-9 uM tRNA transcripts and 3-300 nM 
aaRSs. 

Phylogenetic analyses 

For phylogenetic analysis of mitochondrial tRNA se- 
quences, tRNAs were predicted in mitochondrial 
genomes of interest with Erpin/RNAweasel (21,22) and 
a tRNA profile specific for fungal mitochondria that 
may be used via our website (http://megasun.bch 
.umontreal.ca/RNAweasel). Sequence columns represent- 
ing the sequence positions in the variable loop between the 
anticodon- and T-loop were removed as the nucleotides in 
this region are extremely diverse in sequence and length, 
and phylogenetic analysis was conducted using either a 
Bayesian [PhyloBayes; (23)] or maximum likelihood 
approach [RaxML with GTR model; (24)]. For inference 
of a species tree, an application developed in-house 
(Mams) was used for automated mitochondrial protein 
alignment, removal of ambiguous regions in the align- 
ment, and concatenation. Briefly, derived Cob, Cox 1-3, 
Atp6, 9 and Nad 1-6 protein sequences are pre-aligned 
with Muscle (25), alignments are iteratively refined with 
HMMalign (S. Eddy, http://hmmer.janelia.org) using E- 
values obtained with respective HMM models as an opti- 
mization criterion. Sequence positions that are not aligned 
with a posterior probability value of 1 are removed from 
the resulting alignment. This dataset with 22 taxa and 
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Figure 1. Nucleotide sequences of S. cerevisiae mitochondrial tRNAs. (A) tRNA| 1 with an 8-nt anticodon loop and a UAG anticodon. 

(B) tRNAl hr with a canonical UGU anticodon. (C) tRNA Hls with a G_[ (circled, a major anti-determinant for ThrRS). The primary sequences 
of tRNA[ Br and tRNA His are 72% identical (shaded). 
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3728 amino acid positions was analyzed using Bayesian 
inference by PhyloBayes that implements the 
CAT + GTR model, known to be among the least sensitive 
to LBA artifacts [(23,26,27), and references therein]. 



RESULTS 



Phylogenetic analysis reveals that tRNA™ r originated 
from tRNA His 

The unsolved question of the yeast mitochondrial CUN 
codon reassignment is the evolutionary origin of tRNA] 111 ' 
with an 8-nt anticodon loop and a UAG anticodon. To 
analyze the tRNA] 111 recruitment, we performed a phylo- 
genetic analysis of all mitochondrial tRNAs of 
S. cerevisiae and related yeast species using Bayesian in- 
ference (Figure 2). Not surprisingly, the 10 organisms with 
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Figure 2. Phylogeny of yeast mitochondrial tRNAs. The phylogenetic 
analysis with PhyloBayes (default model parameters) contained all 
tRNA sequences from the species shown in Figure 3. Only the 
section of the tRNA phylogeny covering the tRNA} hr and tRNA Hls 
clusters are shown (marked red and blue, respectively), confirming 
monophyly of tRNAf ' and a sister group relationship to tRNA Hls . 
The posterior probability support for the two tRNA groups is 1.0 
(note that phylogenetic analysis with tRNA sequences depends on 
only few informative nucleotide positions, which does not allow to 
resolve the branching order within these groups). Removal of the anti- 
codon sequence positions from the dataset did not change clustering 
into the two tRNA groups, nor did an analysis of the same datasets 
with maximum likelihood and the GTR model. 



tRNA]" hr are closely related; yet this tRNA 1 " hr cluster is 
related to tRNA Hls (72% sequence identity between the 
two respective S. cerevisiae tRNAs; Figure 1). In contrast, 
tRNAj and tRNAj hr are definitively distant in the phyl- 
ogeny (Figures 2, Supplementary Figure SI) and the re- 
spective 5. cerevisiae tRNAs share only 52% sequence 
identity. As the phylogenetic signal in tRNA sequences 
is limited by the small number of informative sequence 
sites, and not a reliable marker in species phylogenies 
because of occasional identity shifts [e.g. (28)], we have 
built a yeast species tree based on mitochondrion-encoded 
proteins to permit mapping of evolutionary changes 
to this tree. The result of this phylogenetic analysis 
(Figure 3) is consistent with a single origin of the 10 
yeast species that have a tRNA] 1 " homolog with an 8-nt 
anticodon loop, a clade emerging close to the divergence 
of Kluyveromyces species and subsequent to Pichia 



canadensis. Together, these results strongly suggest that 
tRNA} hr derived from mitochondrial tRNA fil: 
ancestral yeast species. 



m an 



Saccharomyces cerevisiae MST1 threonylates both 
tRNA™ r and tRNA^ 

The mitochondria of all yeast species containing tRNA] 111 
have lost the tRNA^ G gene, leaving tRNA, ' the only 
tRNA capable of reading CUN codons. Whereas 
5. cerevisiae mitochondrial LeuRS (ScmtLeuRS) effi- 
ciently attaches Leu to tRNA^ A (29,30), it could not 
recognize tRNA] 111 ' (Figure 4A), confirming that the 
S. cerevisiae mitochondrial CUN codons are assigned to 
Thr instead of Leu. It is known that the WT S. cerevisiae 
strain attaches Thr to tRNA]" hr , and an MST1 deletion 
strain fails to threonylate tRNA] 111 , resulting in a 
respiration-deficient phenotype (16). This suggests that 
MST1 is a putative mitochondrial ThrRS. Compared 
with bacterial ThrRSs, yeast MST1 enzymes lack an 
N-terminal editing domain, but share homologous cata- 
lytic and tRNA anticodon binding domains. To test 
aminoacylation of tRNA] 111 directly, the 5. cerevisiae 
MST1 gene was cloned into pET28a for protein 
overexpression in E. coli. N-terminal His-tagged MST1 
was purified to homogeneity and tested in aminoacylation 
reactions in the presence of [ 3 H] Thr and in vitro 
transcribed S. cerevisiae tRNAs. Consistent with previous 
in vivo results (16), ScMSTl was able to charge tRNA[ hr 
with Thr in vitro (Figure 4B). MST1 could also 
threonylate tRNAj hr , which was unexpected as it was pre- 
viously thought that a second mitochondrial ThrRS was 
responsible for tRNAj' 11 aminoacylation. Steady-state 
kinetic experiments revealed that ScMSTl recognized 
tRNA]" hr and tRNA™ 1 ' with high affinity, with K m 
values of 0.29 and 0.44 uM, respectively (Table 1), suggest- 
ing that tRNA modifications (17) are not critical for 
MST1 recognition. These data, together with previous 
in vivo results, establishes unequivocally that tRNA] 1 " is 
indeed threonylated by MST1. While we favor that MST1 
threonylates tRNA?' 11 ' in vivo, we do not exclude the pos- 
sibility that a second ThrRS activity is present in yeast 
mitochondria. 
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Figure 3. Pliylogeny of yeast species based on concatenated mtDNA-encoded protein sequences. The phylogenetic analysis with PhyloBayes and the 
CAT model is based on 13 mtDNA encoded proteins. All divergence points are supported by posterior probability values of 1.0, except where 
indicated. The red arrow points to the concomitant loss of all seven nad genes and the start of mitochondrial codon reassignments, including AUA 
methionine, CUN threonine. Species shown in black possess mitochondrial tRNA[^ G (but not tRNA^ hr ), and in these organisms CUN codons are 
assigned to Leu. The yeast species marked red, such as K. thermotolerans, have lost mitochondrial tRNAj^ G and obtained tRNA™ 1 that decodes 
CUN codons as Thr. K. lactis is marked magenta as it has no CUN codons and no corresponding tRNA with a UAG anticodon. A. gossypii 
(marked blue) contains a tRNA species with a UAG anticodon, yet its identity of this tRNA and the amino acid reading CUN codons in A. gossypii 
remain obscure (to be discussed by BFL in Organelle Genetics: Evolution of Organelle Genomes and Gene Expression, Springer 2011). 
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Figure 4. Aminoacylation by S. cerevisiae 
(B) Threonylation of tRNA, 

and tRNA His variants (3 uM) by ScmtHisRS (0.3 uM). 



MST1 and HisRS. (A) Leucylation of tRNA™ 1 " (3uM) and tRNA^A (3 uM) by ScmtLeuRS (0.3 uM). 



(3 uM), tRNAT hr (3 uM) and tRNA His variants (3 uM) by ScMSTI (0.3 uM). (C) Histidylation of tRNAf hr (3 uM) 



Loss of G_i converts S. cerevisiae tRNA ls to a substrate 
for S. cerevisiae MST1 

To provide experimental evidence for the evolution of 
tRNA}' 11 ', we examined recognition of tRNA Hls variants 
by MST1. Except in a few a-proteobacteria [(31), and ref- 
erences therein], all known tRNA Hls species contain a G 
at position —1, which is a critical identity element for 
histidyl-tRNA synthetase (HisRS) [(32,33) and references 
therein]. Sequence alignments of tRNA]" hr and mito- 
chondrial tRNA Hls genes revealed that G_i addition in 
tRNA Hls comprises one major difference between the 
two tRNA species (Figure 5). In vitro, ScMSTI failed to 



charge the WT SctRNA His transcript with Thr, but 
deleting G_i stimulated threonylation of tRNA Hls by 
ScMSTl (Figure 4B). Steady-state kinetic data showed 
that AG_i tRNA His gained 4% threonylation activity of 
the WT tRNA™ 1 ' (Table 1). Compared with tRNA}' 11 ', 
AG_| tRNA His displayed 10-fold lower k cat and 3-fold 
higher K m values for threonylation by ScMSTI . In 
addition to G_i, other major differences between 
tRNA His and tRNA™ r include an A insertion in 
tRNA^ hr at position 35, and the discriminator base at 
position 73 (Figure 5). Changing C73 to A did not 
improve the threonylation efficiency of AG_i tRNA Hls , 
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Table 1. Threonylation by SfMSTl 

Anticodon loop k CM (min~') K m (\iM) k CM jK m (uMmin -1 ) Relative k all jK, 



tRNA™ r 












WT 


UUUUAGGU 


2.8 ± 0.4 


0.29 ± 0.09 


10.3 ± 2.3 


100 


G_, 


UUUUAGGU 


0.094 ± 0.032 


1.5 ± 0.6 


0.065 ± 0.006 


0.63 


tRNAj hr 












WT 


UUUGUAA 


2.3 ± 0.1 


0.44 ± 0.04 


5.4 ± 0.6 


52 


tRNA His 












WT 


UUGUGGU 


ND 


ND 


ND 


ND 


InsA35 


UUGUAGGU 


ND 


ND 


ND 


ND 


InsA35/C73A 


UUGUAGGU 


ND 


ND 


ND 


ND 


AG_| 


UUGUGGU 


0.27 ± 0.10 


0.74 ± 0.41 


0.39 ± 0.06 


3.8 


AG_,/InsA35 


UUGUAGGU 


0.34 ± 0.07 


0.20 ± 0.02 


1.7 ± 0.5 


17 


AG_,/C73A 


UUGUGGU 


0.23 ± 0.01 


0.85 ± 0.32 


0.29 ± 0.09 


2.8 


AG_,/C73A/InsA35 


UUGUAGGU 


1.2 ± 0.4 


0.57 ± 0.26 


2.4 ± 1.0 


23 



ND, not determined due to low activity. 
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tgaatatatttcaat-ggt-ag-aaaagacgcttgt 
:taaatatatttcaat-ggttagcaaaatacgcttgt 
3tgaatatatttcaat-ggt-ag-aaagtatgcttgt 
tgaatatatttcaat-ggt-ag-aaaatacgcttgt 
taaatatatttcaat-ggt-ag-aaaatatgcttgt 
tgaatatattttaat-ggt-aa-aaagtacgcttgt 
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tgaatatatttcaat-ggt-ag-aaaatacgcttgt 
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taaatatatttcaat-ggt-ag-aaaaagtacttgt 
taaatataatttaacaggt-aa-aatgtatgttttt a 
tagatataatttaatcggt-aa-aatgtatgttttt a 
taaatatagtttaat-ggt-ag-aatatatgttttta 
taaatataatttaat-ggt-aa-aatgtatgttttt a 
taaatataatttaatcggttaa-aatgtatgttttta 
taaatataatttaat-ggttaa-aatatatgttttt a 
taaatataatttaat-ggt-aa-aatgtatgttttt a 
taaatataatttaat-ggt-aa-aatgtatgttttt a 
taaatataatttaataggt-aa-aatgtatgtttct a 
taaatataatttaat-ggt-aa-aatatatgttttt a 
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'.gtgcgtttaatctaagttcaattcttagtattcace 
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Figure 5. Sequence alignment of mitochondrial tRNA™ r and tRNA Hls . Three major differences between tRNA, lrir and tRNA"'" sequences are 
indicated by boxes. Cc, Candida castellii; Cg, Candida glabrata; Kt, Khiyveromyces thermotolerans; Nb, Nakaseomyces bacillisporus; Nd, 
Nakaseomyces delphensis; Sea, Saccharomyces castellii; Sc, Saccharomyces cerevisiae; Sp, Saccharomyces pastorianus; Ss, Saccharomyces servazzii; 
Vp, Vanderwaltozyma polyspora. 
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but inserting A35 in the anticodon loop of AG_i tRNA Hls 
further increased its threonylation activity 5-fold 
(Table 1). In the presence of G_i, A35 insertion did not 
allow threonylation of tRNA Hls . These results suggest that 
G_i is a major anti-determinant in tRNA Hls for MST1. 
In line with this notion, addition of G_i to tRNA| hr 
reduced its threonylation activity 150-fold (Table 1). 

G_i addition allows tRNA[ hr to be histidylated by 
S. cerevisiae HisRS 

ScMSTl efficiently aminoacylates mitochondrial tRNA Hls 
variants lacking G_ i . We next investigated the recognition 
of tRNA His and tRNA} hr variants by mitochondrial 
HisRS. In S. cerevisiae, a single nuclear gene HTS1 
containing two open reading frames encodes both the 
cytoplasmic and mitochondrial forms of HisRS (34). 
The mitochondrial HisRS (mtHisRS) harbors an extra 
mitochondria-targeting sequence, yet the predicted 



mature form of the mitochondrial enzyme is identical to 
the cytoplasmic HisRS. Previous studies revealed that 
ScHisRS recognizes both the G_j and the anticodon of 
the cytoplasmic tRNA Hls (35,36). We overexpressed 
mature ScmtHisRS in E. coli and purified the enzyme to 
homogeneity. In vitro aminoacylation reactions showed 
that ScmtHisRS efficiently charged His to the WT mito- 
chondrial tRNA His (Figure 4C, Table 2), with a K m value 
comparable to that of the cytoplasmic tRNA Hls (35). 
Surprisingly, neither deleting G_i nor inserting A3 5 in 
the anticodon loop of mitochondrial tRNA Hls significant- 
ly affected the aminoacylation activity. However, when 
both changes were introduced into mitochondrial 
tRNA Hls , the histidylation activity decreased by 80-fold. 
The discriminator base also appeared to be important for 
recognition by ScmtHisRS. These results indicate that in 
the context of S. cerevisiae mitochondrial tRNA Hls , G_i is 
dispensable for HisRS recognition and the anticodon can 
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alternatively serve as the major identity element. 
ScmtHisRS was unable to aminoacylate WT tRNAT hr 
(Figure AC). Interestingly, addition of G_! to tRNA? ir 
restored 7% histidylation activity of the WT tRNA his 
(Figure 4C, Table 2), further strengthening our 
argument that the mitochondrial tRNA His and tRNA[ hr 
in yeast are closely related during evolution. 

MST1 has coevolved with tRNA[ hr to establish CUN 
codon reassignment in yeast mitochondria 

The biochemical and phylogenetic evidence above 
suggests that tRNA™ 1 has evolved from mitochondrial 
tRNA Hls , at a time point close to the divergence of 
Kluyveromyces species. To understand whether tRNA™ 1 
could be recognized by MST1 enzymes from other fungal 
species, we overexpressed and purified MSTls from 
C. albicans and S. pombe [a non-hemiascomycete 'fission 
yeast' belonging to Taphrinomycotina (37)], and tested 
them for threonylation of tRNA} hr . CaMSTl and 
SpMSTl share 49 and 43% sequence identity with 
ScMSTl, respectively, and both enzymes were able to 
aminoacylate S. cerevisiae mitochondrial tRNAj hl 
(Figure 6A). However, neither CaMSTl nor SpMSTl 
recognized tRNA™ 1 ' or the tRNA Hls variants tested 
above (Figure 6). These results strongly suggest that 
MST1 specifically evolved to recognize tRNA™ 1 in a 
group of yeasts comprising S. cerevisiae. Therefore, 



CUN codon reassignment was completed following the 
coevolution of MST1 and tRNA™', which established 
specific protein-tRNA interactions. 

DISCUSSION 

Reassignment of CUN codons from Leu to Thr in yeast 
mitochondria 

CUN codons in the mitochondria of Saccharomyces, 
Nakaseomyces and Vanderwaltozyma have been previously 
assigned to Thr (5,11,12). This recoding event is ac- 
companied with the loss of tRNAj}™ G and the appearance 
of tRNA™ 1 with an unmodified UAG anticodon (17). 
Analysis of a mass spectrometry database (PeptideAtlas) 
confirms that CUU and CUA codons indeed encode Thr in 
S. cerevisiae mitochondria. CUG and CUC codons are 
rarely used in yeast mitochondria, and the nature of the 
amino acid translated by such codons has not been 
verified experimentally. Given that the unmodified U at 
the first anticodon position is able to recognize all 4nt 
(6,38), and no tRNA bearing GAG or CAG anticodons 
have been shown to be encoded by or imported into yeast 
mitochondria (39), it is plausible that CUG and CUC are 
also decoded by tRNA™ r as Thr. tRNA™' possesses 
an enlarged anticodon loop with 8nt (UUUUAGGU), 
whereas a canonical tRNA anticodon loop consists 



Table 2. Histidylation by ScmtHisRS 

Anticodon loop k cai (min~') K m (uM) kat/Km ( U M min~') Relative k a JK, 



tRNA7 hr 












WT 


UUUUAGGU 


ND 


ND 


ND 


ND 


G_, 


UUUUAGGU 


1.8 ± 0.6 


0.43 ± 0.20 


4.3 ± 0.6 


7.2 


tRNA H,s 












WT 


UUGUGGU 


25 ± 7 


0.43 ± 0.18 


60 ± 10 


100 


InsA35 


UUGUAGGU 


27 ± 8 


0.60 ±0.13 


45 ± 4 


74 


InsA35/C73A 


UUGUAGGU 


38 ± 3 


0.28 ± 0.01 


130 ± 16 


220 


AG_i 


UUGUGGU 


14 ± 3 


0.40 ± 0.22 


45 ± 25 


75 


AG_,/InsA35 


UUGUAGGU 


1.8 ± 0.1 


2.4 ± 0.8 


0.78 ± 0.22 


1.3 


AG_,/C73A 


UUGUGGU 


1.2 ± 0.2 


0.54 ± 0.28 


2.6 ± 1.0 


4.3 


AG_i/C73A/InsA35 


UUGUAGGU 


0.75 ± 0.07 


2.2 ± 0.5 


0.35 ± 0.06 


0.58 



ND, not determined due to low activity. 




CaMSTl /SpMSTl 



-□-CaMSTl + tRNA 1 Thr 
-O- CaMSTl + tRNAJ hr 
SpMSTl + tRNAJ 1 * 
SpMSTl + tRNA™' 



B 

1.6 - 

1.2 • 

3.8 - 

).4 • 
0 

10 0~ 



CaMSTI 



-□-WT tRNA His 
-O-AG. , tRNA His 
-&-AG ,/C73A tRNA His 
-^-AG ,/C73A/lnsA35 tRNA His 



c 
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■\ 0.8 
0.4 



SpMSTl 



4 6 
Time (min) 



8 



-□-WT tRNA His 
-O-AG., tRNA"» 
-&-AG.,/C73A tRNA His 
-V-AG.,/C73A/lnsA35 tRNA H,s 



Figure 6. Threonylation of tRNAf hr and tRNA His variants (3 uM) by C. albicans and S. pombe MST1 (0.3 uM). (A) CaMSTl and SpMSTl 
threonylate tRNAj hr but not tRNA[ hr . (B and C) CaMSTl and SpMSTl are unable to threonylate tRNA His variants. 
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of 7nt with the triplet anticodon centered at positions 
34-36. In addition to tRNA]" hr , several natural tRNAs 
have been found to contain eight bases in the anticodon 
loop, yet such tRNAs decode quadruplet codons instead of 
triplet codons (40). A question that arises is whether 
tRNA]" hr can also decode quadruplet codons. In yeast 
mitochondria, several protein-encoding genes, such as 
coxl, harbor in-frame CUAA sequences. Decoding of 
CUAA by a quadruplet anticodon could cause detrimental 
frameshift translation of these critical mitochondrial 
proteins. It is thus reasonable to think that tRNA]" hr 
decodes all four CUN triplet codons but not quadruplet 
codons. 

Evolutionary origin of tRNA™ r 

Our analyses reveal that tRNA^ hr was recruited from 
mitochondrial tRNA Hls . It has been long considered that 
tRNAs accepting the same amino acid (isoacceptors) 
evolved by gene duplication from the same common 
ancestor. However, studies in bacteria demonstrate that 
several tRNA genes may derive from different amino 
acid accepting groups (alloacceptors) (41). Naturally- 
occurring recruitment of alloacceptor tRNA genes was 
later reported in sponge mitochondria based on phylogen- 
etic studies (28). Our work combines both biochemical 
and phylogenetic approaches to provide the first clear 
evidence that alloacceptor tRNA gene recruitment is 
directly responsible for a recent recoding event, suggesting 
that such a mechanism might have played an important 
role in the establishment of the genetic code during 
evolution. 

Phylogenetic analysis suggests that tRNA| hr likely oc- 
curred after the split between Candida and Saccharomyces, 
subsequent to the divergence of P. canadensis. Interestingly, 
the C. albicans and C. parapsilosis mitochondrial genomes 
each contains two copies of the tRNA Hls gene but no 
tRNA™ 1 , and the CUN codons in these mitochondria 
remain assigned to Leu. In Kluvyeromyces lactis, the mito- 
chondrial genome lacks both CUN codons and any tRNA 
reading such codons (7). Previous studies also suggest that 
current CUN codons decoded as Thr did not originate from 
ancestral Leu codons, but instead are derived from codons 



for other amino acids (18). Although we do not exclude 
the possibility that the original CUN codons could be 
ambiguously assigned to both Leu and Thr in certain inter- 
mediate yeast species, we favor that the emergence of 
tRNA'[ hr is a late event following the loss of original CUN 
codons, and such codons have reemerged in present-day 
yeast mitochondria. We propose that CUN codon reassign- 
ment initiated with duplication of the tRNA Hls gene, 
followed by the loss of CUN codons and tRNA[}™ G 
(Figure 7). Next, one copy of the tRNA Hls gene evolved to 
the present-day tRNA]" , and MST1 coevolved with 
tRNA]" 111 to form a cognate aaRS/tRNA pair. Finally, the 
CUN codons reappeared from various other codons and 
were reassigned to Thr. Our biochemical studies suggest 
that HisRS could still recognize the intermediate tRNA 
during the evolution from tRNA His to tRNA[ hr (Table 2). 
However, the mitochondrial proteome remained largely un- 
affected as the CUN codons were not present in the mito- 
chondrial genome during this period. 

Several mechanisms have been proposed to explain 
the codon reassignment processes, including: (i) codon 
capture mechanism (42), hypothesizing that a specific 
tRNA and corresponding codons completely disappear 
from a genome before a novel tRNA evolves to read 
such codons; and (ii) ambiguous intermediate mechan- 
ism (43), postulating that a tRNA mutant is recognized 
by more than one aminoacyl-tRNA synthetases, and a 
codon may be assigned to multiple amino acids in 
an intermediate state. The latter mechanism is supported 
by the discovery that in several extant Candida species, 
the CUG codon is ambiguously decoded by both Leu 
and Ser due to the presence of a tRNA CAG charged by 
both leucyl- and seryl-tRNA synthetases (12,44). Our 
work suggests that CUN codons and tRNA{}™ G are lost 
in yeast mitochondria prior to the emergence of tRNA| hr , 
which provides a paradigm for the codon capture mech- 
anism and lends support to the evolving genetic code 
theory. 
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Supplementary Data are available at NAR Online. 



LeuRS HisRS MST1 LeuRS HisRS MST1 HisRS MST1 




Figure 7. Proposed model for CUN codon reassignment in yeast mitochondria. (1) tRNA Hls duplicated in an ancestral yeast species while CUN 
codons remain assigned to Leu. (2) CUN codons were changed to UUA or UUG decoded by tRNA^' A , and tRNAj}™ G was lost. (3) tRNA lb 
evolved to tRNA™ 1 , and CUN codons reemerged from various codons. 
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